The nonsingular brane solutions via the Darboux transformation 



o 
o 

(N 

Oh: 
< 



A.V. Yuro\E and V.A. Yuro\|l| 

The Theoretical Physics Department, 
Kaliningrad State University, 
A. Nevskogo str., 14, 236041, Russia. 

(Dated: February 1, 2008) 

We consider the Darboux transformation as a method of construction of exact nonsingular so- 
lutions describing the three-dimensional brane that interacts with five-dimensional gravity and the 
bulk scalar field. To make it work, the five-dimensional Einstein's equations and the Israel's con- 
ditions are being reduced to the Schrodinger equation with the jump-like potential and the wave 
functions sewing conditions in jump point correspondingly. We show further that it is always pos- 
sible to choose the functions in Crum's determinants in such way, that the five-dimensional Ricci 
scalar R will always be finite both on brane and in bulk. The new exact solutions being the gen- 
eralizations of the model with the odd superpotential are presented. Described formalism is also 
appliable to the cases of more realistic branes with cosmological expansion. As an example, via the 
usage of the simple orbifold model (5*1 /Z2) and one-time Darboux transformation we construct the 
models where the cosmological constant on the visible brane is exponentially small. 

PACS numbers: 11.25.Mj, 11.27.-(-d 
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I. INTRODUCTION 

In article T| , there have been considered an interesting 
model, which describes the set of parallel 3-branes, all of 
them being imbedded in five-dimensional bulk, filled with 
the gravitation and the scalar field. The action of set is 
expressed as: 



S = J d^xdy 




(1) 



where {x^,y) and < /i < 3 are the five-dimensional 
coordinates; brane number b is located at y = ?/(,; gah is 
the 5-D metric, and g^^^ is the metric, induced on the 6th 
brane. (Note that it is written in special system of units, 
where the gravitational constant n = \ and speed of light 
c — 1.) The tension on the &th brane is named (Tb, and 
the potential V{4)) is considered to be a function of the 
bulk scalar field (j) = 4'{y)- The field equations follows 
from the Q by the procedure of variation; their solution 
is assumed to have a form (Q): 



dy' + A{y) (- 



-dt^ + e^"*S,kdx'dx'') , 



(2) 



where A{y) is a sought warp factor, and iJ=const is a 
Hubble parameter, li H = on the brane, the brane will 
be the stationary one, and if 7? > 0, then the brane will 
be expanding in the de Sitter regime 

The (P) leads to system of differential equations, that, 
in turn, can be rewritten in "supersymmetric" style by 
introduction of the superpotential W{(j)), being defined 
by the (j) = W'{(j))/2 relation, where the dot denotes the 



derivative in y, and the prime - the derivative in (j). In 
case of ff = there exists following equation 



V{(f>) = M/'V8- VFV6. 



(3) 



Upon usage of this equation, the authors of H have 
suggested a simple way of obtaining of exact solutions of 
the field equations (see also Q, 3 for further details). 
Their approach can be expressed in following steps: 

(i) Choosing convenient V{(j)) and W{(j)), making sure 
that the equation Q with these new values still holds 
true; 

(ii) Integration of equation = W'{(t))/2, in order to re- 
ceive the relation </> = 4>{y)', 

(iii) When H = there exists the relation A{y)/A{y) — 
—W{(j){y))/3. Hence, substitution of the obtained func- 
tion for (f>{y) in explicit equation for W{(l)) and further 
integration results in the sought warp factor A{y). 

From this point, the famous solutions RS (Randall- 
Sundrum) turns out to be just a case W{(/)) = zLa^ = 
const. Authors U have also considered the models with 
even iW{(t)) = ± (a^ - bcj)^)), odd {W{(j)) = ±2b^(t>) and 
exponential {]¥{(/)) = 2ae~^'^) potentials (the even su- 
perpotentials had been previously discussed in [5, and 
the odd ones had been introduced in 0). 

The major advantage of the described method lies in 
it's simplicity. For example, it allows an easy construc- 
tion of nonsingular solutions (see also 0). Lets consider 
the case of one brane, localized in y = j/o- From the equa- 
tion (j) = W'{(f))/2 it is possible to express y as a function 
of the field variable 6: 



4>o 



W'{(l))' 



(4) 
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where (jjo — 4>{yo) is a value of a scalar field on brane. 
It follows from the Q that the very existence of zeros 



2 



of the function W'{4>) {4> = (j)^ = const, | |< oo) 
is leading to the divergences of the integral. In another 
words, (j) ^ (fy^ when ?; — > ±cx). What it means, is that if 
we choose the superpotential as having at least two zeros 
at both (p > <j>o and < 0o where W' is vanishing fast 
enough (for example W'{(t)) ^\ 4> ~ with k > 1, 

and <j) ~^ 4'^)j then the field remains finite throughout 
all the bulk space. On the other hand, the case of H ~ 
allows to express the Ricci scalar via the superpotential 
and its derivative in the form: 



Ricj,)^2iW'mV3-5W\<P)/9, 



(5) 



thus, showing the finiteness of the scalar curvature for 
the aforementioned nonsingular superpotential for every 
y. The particular examples of the nonsingular solutions 
for even, odd and exponential superpotentials can all be 
found in 0. It also contains the generalization of this 
technique for the non-stationary branes {H ^ 0) and 
a brief discussion touching the possible usages of the 
method in solving such a well known problems as the 
problem of cosmological constant's smallness and the hi- 
erarchy problem. 

However, in spite of all the strong points, this method 
does have its apparent weaknesses, the primal one be- 
ing its " phenomenological" nature - the fact, that each 
physically consistent (i.e. leading to the reasonable (f>{y) 
and A{y)) superpotential W{(t)) can only be obtained by 
"guessing". Although all simple examples, considered in 
[l| were not too hard "to guess" and they do lead to 
sensible cosmological models, it would not necessarily be 
that easy, whenever we try the more difhcult (and hence, 
more realistic) W{(f>). These reasonings lead us to prin- 
cipal question: does there exist such procedure of (^j's 
nonsingular solutions construction, which is both math- 
ematically simple and regular? In this work we give the 
positive answer for this question and introduce the new 
method, which satisfies both of the conditions. 

The general idea of our approach can be formulated 
as follows: the supersymmetric link between V{(f)) and 
14^(0) is a typical one for the supersymmetric quan- 
tum mechanics (SSQM). In turn, SSQM is realizable 
through the usage of Darboux transformation (DT) for 
the Schrodinger equation 0, and that allows us to 
construct the integrable nonsingular potentials - provid- 
ing that all the prop functions of Crum's determinants 
do have the intermittent zeros 0, So, it seems to be 
only natural trying to write the DT for the field equations 
directly, and then using them as a source of a new exact 
solutions. An existence of quite simple algebraic con- 
nections between the Ricci scalar and the superpotential 
(O notes that construction of nonsingular potentials in 
SSQM via the right DT application will also aid in devel- 
opment of the systematical procedure of construction of 
the nonsingular solutions (i? ^ oo) in the brane theory, 
which, by that means, will turn out to be a regular one - 
as distinct from the " phenomenological approach" of . 
In particular, the DT method, if applied to just a single 
exactly solvable potential, grants us the infinite number 



of an exact nonsingular solutions, with the mentioned po- 
tential being a " fuse" in all those solutions construction. 
For example, in this work we'll take the RS for the role 
of such a fuse solution. Another interesting way of the 
DT usage is a construction of the chains of discrete sym- 
metries and the further inquiries of their closures. Such 
a technique allows getting the new exact solutions, which 
are manifested through the Painleve transcendents and 
the higher extensions. 

The article has the following structure. Next para- 
graph is dedicated to the Darboux transformations and 
to their compatibility with the Israels jump conditions 
at the brane j9|- In third paragraph we'll present some 
new exact solutions, received by one-time Darboux dress- 
ing of the RS solutions. Double DT, the importance of 
the shape-invariance and the dressing chains are all ex- 
amined in the fourth paragraph. In particular, it will 
be shown, that even, odd and exponential superpoten- 
tials, discussed in work ^ are just the examples of the 
shape-invariants with regard to DT potentials. In the 
same section we'll show the exact solutions, that serve 
as the generalizations of the model with the odd super- 
potential (here we include the solutions expressed via 
the fourth Painleve transcendent). In paragraph 5 we'll 
describe the construction of an exact solutions for the 
non-stationary branes, made with the usage of DT and 
will discuss the nascent difficulties. Here we'll consider 
both noncompact and compact spacetimes. In the last 
case, using the simple orbifold model {S1/Z2) and one- 
time Darboux transformation we'll construct the models 
with two branes where the cosmological constant on the 
one of them (visible brane) will be exponentially small. 
Also, we'll outline the new approach (called Green func- 
tion method) which allows to construct solutions both 
in stationary and non-stationary cases. Finally, the six 
paragraph is an overall conclusion. 



II. DARBOUX TRANSFORMATIONS AND 
JUMP CONDITIONS 

li H = and u = A/A, then from the (QJ and (|2J) we 
get the following set of equations ilj: 

~^-^^cFb{(l})Siy-yb), w2 = y-^F(0). 



Third equation has the form 

91/(0) , ^ dab{(t>) 



(6) 



(j) + 2u4> : 



and is not independent, because it can be received from 
the set (01. In case of the single brane, located at y = 0, 
we get a{(j){y))S{y) = (Jo5{y), where ctq = o-(0(O)) (sup- 
posing 0(0) 7^ 0). It is easy to check, that the function 
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A'^{y) satisfies the Schrodinger equation 

(^ + ^my))+/3%))^'(2/) = 0, (7) 

where /3 — 4cto/3. Now lets define the spectral problem 

i^^iv{y)-pS{y)-X)^, (8) 

where ip = ipiy'i A) and v{y) is a certain potential which 
will be further refereed to as a quantum-mechanical po- 
tential, or QM-potential. Lets assume, that it is possible 
to solve the equation ((SJ formally (i.e., obtaining among 
others some "non-physical" solutions which do not be- 
long to L^) for some QM-potential v{y) and all A G M. 
Lets then suppose that for some values of the spectral 
parameter A (say, A = Aq) the solutions of the lO will be 
strictly positive: V'(2/i Ao) > 0. In this case we can con- 
clude that function A(y; Aq) = ipiVi Ao) is the solution 
of the with potential 

my)) = ^(Ao-t'(y)). (9) 

Therefore, the equation ((SJ can, in fact, be considered as 
the producer of exact solutions of the ||7J). Substituting 
A{y) and V((j){y)), obtained from the ||SJ|, in ^ we'll 
find out 4> = (j){y). In this point there can be raised the 
question: what solutions can be named the physically 
sensible one? More precisely, what can be said about 
the solutions with (jp < 0? 

Nota Bene 1. Recently the great deal of interest 
has been drawn to those very specifical cosniological 
models, which adrnits the violation of the dominant- 
energy condition |jlo|. Such models often presuppose 
that the "dark energy" has so-called phantom compo- 
nent which should be described by some scalar 
field with the negative kinetic term. Quantum theory of 
those fields has been considered in [13| (this work also 
has an almost complete bibliography of a subject). As a 
general result of this investigations we can point out the 
following indication: apparently such fields fails to be 
consistently considered as a fundamental objects. But 
nevertheless, we can't totally exclude the probability 
that the Lagrangians with the negative kinetic terms 
will somehow appear in a role of effective models. 
For example, similar terms appear in some models of 
supergravityng and in the gravity theories with highest 
derivatives Finally, the "phantom energy" in the 
brane theory was considered in [l5j |. However, despite 
of all the above said, we will not further consider the 
phantom component. In other words, we shall suppose 
that (p^ > during the rest of this paper. 

For specificity, we'll further assume the following con- 
dition for the potential v{y): 

/ + 00 
dy (1 + \y\) \v{y) < oo. 
-oo 

(10) 



It is also worth mentioning, that for our further inves- 
tigation the validity of (|10|l is not really necessary. We 
just found it to be convenient to consider the potentials, 
that arc confined within such bounds. 

So, lets assume the parameters tp{y)i ''^(y)j A to be al- 
ready determined. We are going to find the functions 
A{y), <j){y) and V{(l)). In order to do this we'll examine 
the region located far from the brane to exterminate the 
(5-function term in the © from further considerations. 
It's obvious, that the quantity A(jj) can only be restored 
if ■0(2/) is nonnegative for all y. The only one wave func- 
tion that both satisfies this condition and belongs to the 
physical spectrum of Schrodinger equation (jS)) , is a wave 
function of ground state with A = A„^^ . We should also 
mention, that there exists another possible way to choose 
the "good" wave function, namely, to consider the case 
A = A_i < A^^^. This satisfies the imposed condition, 
since the ^ indeed admits the always positive solutions 
ip{y; A_i). These solutions obviously do not belong to the 
L2 space and definitely do n ot descri be the bound states, 
but the warp factor A — \J'ip{y] AJ will nevertheless be 
real. However, these solutions result in singularity at in- 
finity: A — > cxD at y ^ cxD, and therefore the ground state 
-0(j/; A^^J is by far the best alternative to V'(2/; A_i) 

It is well known, that DT has a remarkable capac- 
ity to enable the engineering of ad hoc potentials with 
arbitrary finite discrete spectrum, thus it can be consis- 
tently considered as a procedure of construction of the 
system ((Sjl's exact solutions. We'll now quickly remind 
the very essence of DT method 0. Let tpi = ^{y; Ai) 
and ip2 = V'(y; A2) be the two solutions of the equation 
©, and let also ^pi > for every y it is defined on. The 
"01 is called the prop function of Darboux transformation. 
In general, the transformation law is: 

(1) _ 02-01 - 01-02 



X; ^ ^(l) = W - 2— r log'01, A2 ^ A2. 

Darboux transformation Hll|) is an isospectral symmetry 
of equation because the dressed function 02^'' is the 
solution of the dressed equation : 

^W = (.W-/3W%)-A2)4'\ 

with a new (dressed) potential w'^^ but with the same 

value of Aj^-* = A2. The transformation law for the prop 
function is given by the relation (jlj): 

V'l-V-r^ = iir(Ci + C2/v?(y)d2/), 

(12) 

with arbitrary constants Ci^2 (and, of course, A^^^ = Ai). 
The most interesting part for us here is a ground state 
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of a dressed Hamiltonian. There exist three possibihties. the same conditions (|13|l with /3^^ 



1) If ipi is a wave function of a ground state, and 
A2 is an energy of a first excited level, then the new 
dressed function -02^^ will be the wave function of a new 
Hamiltonian's ground state All other values of the 
discrete spectrum of this Hamiltonian are received by 
simple obliteration of level Ai from the initial spectrum. 

2) The prop function V'l = V'i(?/;'^i) does not belong 
to an L2, but it doesn't has any zeros either, and 



^1 



=+p \y\ 



when y ±00. In the case, after the DT 



we'll get the new Hamiltonian, whose discrete spectrum 
is just an old one except for the new state A = Ai added, 
thereby being the new ground state. The corresponding 
wave function ^[^^ (I12II with C2 = becomes the ground 
state's wave function jlTI] , 

3) The aforementioned prop function can in fact be rep- 
resented as a linear superposition of two positive, linearly 
independent solutions V'l^^ with the following character- 



istics: i/'i^' 



when y =pcx3 and tpi 



(±) 



=+p \y 



y — > ±00. Taking either V'l^^ or ip[ as a prop function 
results in creation of a new Hamiltonian, whose discrete 
spectrum is totally coincides with the initial one .17.]. 
Our article is dedicated in general to the second case. 
As a next step of our consideration we should exam- 
ine the jump conditions, existing on the brane. These 
conditions can be written as 

ip,+io) = ^-,-(0) = MO), V'^+(o)-^.-(o) = -PMo)- 

(13) 

where {y) and ipi- (y) are the solutions of Schrodinger 
equation JSJ, both with the same potential and eigenval- 
ues, but valid only for positive or negative y correspond- 
ingly: 



'4'i± = {v{y) - f35{y) - Xi) ip., 



(14) 



It is easily verified, that (|13(l are exactly the same jump 
conditions as in before-cited article 0: 



K+(0)-j._(0) = i( 



V.+ (0)-i/.-(0) _ 
2V>(0) 



V+(o) ■0_(o) 

2 3 ' 



being, in essence, the Israel conditions 0, 0. The 
next theorem proves, that DT keeps the jump conditions 
(|13|l unchanged: 

Theorem. If iP,{y) = {tP,+ {y > 0),2P,^{y < 0)} 
with i — 1,2 are solutions of (|14|l . satisfying the H13|) 
and if tpi — {"01+ j "01- } is a prop function, then the 

transformed function -02^^ = |'02+' ^'2- 1 (where the 
prop functions of ip^^ are ip-^_^ correspondingly) satisfies 



4+Ho) = 0^-^(0)^ (0), 



(15) 



The proving of the theorem is conducted by the direct 
calculation. We'll just cite two formulas, which turned 
to be quite useful in proving: 



fli = w - Al - Ti±, 



/(I) 



where ti± = ipi±/4'-^±. An obvious result of this the- 
orem is the conclusion, that n-times repeated DT leads 
to the dressed functions, satisfying the (|13|l with tension 
= (— 1)"/3, on the assumption that n prop functions- 
solutions {V-'ii} (* = 1,--,'T-) of the initial Schrodinger 
equation (|14|) and one transformed solution 0ri+i± = 
with the same potential, but different A are all satisfy the 
condition H13() . Such n-times dressing is usually conve- 
nient to introduce via the Crum's formulas Il9l: 



when '^±\y) 



„(") 



where 



^-2 



-1] 



^-2 



-2] 
-l± 



dy'^ 



V'2± 



log A W (2/), 
(16) 



-2] 
-2] 



^-1] 
-i± 



-2] 



and = d'^^'/dy^. These formulas are highly useful for 
studying. In particular, the solution, written in this form, 
taken along with (|15|) . allows us to extend the conditions, 
introduced by Q: 



in). 



n+1 2q"o 
3 ' 



where Uj?-* = yi^"-* 



± /A'l'' and 

(Jo is a tension on the initial 
brane. Moreover, using the H16() in the procedure of levels 
addition one can show, that the wave function of ground 



state for potential Wj?^ is defined by the formula: 



(") - 



(17) 
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Wave function of a ground state is nowhere turns to zero, 
and this statement is similar to condition 7^ 0. Using 
the n-times dressed scale factor, we can calculate the 
corresponding Ricci scalar: 



(") /l(n) 



R 



(n) 



A. 



(n) 



(18) 



where the ± signs are referring to y > (over the brane) 
and y < Q (under the brane). Since the denominator 
of this expression is always non-zero, it follows right 
from the 117|l and p8|l that Ricci scalar remains finite 
not just on brane, but also in a bulk. If, for example, 
the behavior of "01"^ is an asymptotical one, i.e. if 

■(/)(") — > e~P 1^1 when y ±00, then iJj?'' — at 
infinity. Thus, in the case of single stationary brane, the 
DT method allows the simple construction of an abun- 
dant set of an exact nonsingular solutions of equation 

Nota Bene 2. As we have seen, the equation ((S)) 
can be extremely effective in producing of the exact so- 
lutions of ©. In the cases of exactly solvable potentials 
v{y), the usage of DT allows the construction of a large 
(in fact, infinite) set of such solutions. The problem 
arises if the considering Schrodinger equation starts out 
from some arbitrary solution of the ©. This solution 
gives us, in turn, the single exact solution of the {T)) 
with v{y) = —8V{(j){y))/3 and A = - but, generally, it 
is not guaranteed, that v{y) = —8V{(j){y))/3 will be one 
of exactly solvable potentials, and hence - not provided 
that we can solve the spectral problem © for all other 
admitted A. This means, that, in order to use the DT 
method, we should in first place carefully choose the 
initial solution of the ©. It is highly remarkable, that 
the famous RS model actually does the trick! 



III. ONE-TIME DRESSING OF THE RS 
BRANE. 

To illustrate the effectiveness of formulas (|12|) 
we'll take the example of the RS brane dressing, which 
case corresponds to a simple quantum potential v = fi"^ = 
const. Here v = (see (llUll '). and the potential of a 
scalar field V — S^^/S. Solution ipi of equation ^ 
with eigenvalue Ai = has the form 



^P,^=a,^e^^y + b,^e-^^y, Vi+=^i(y>0), 

^p,_ = My < 0), 



(19) 



where a^j^, and b-^^ are real positive constants. Usage of 
(|13ll results in following equations: 



a, = a. 



2^ 



b,_ ^ b,^ 



+&1+) 



(20) 



It is obvious from l|20() that a^^ -|- 6^^ = + b-^_ , there- 
fore it is convenient to take advantage of the parameter- 
ization: 



2-2 I 

= r sm , o^^ 



r^cos^a^. (21) 



(|21|l and H20() allows us to get the compatibility condi- 
tions: 



.2 /? 2 

sm a. < — < cos a, , 
+ 2^ + 



(22) 



where \/3\ < 2fi necessarily. Thus, in order to get the 
new exact solutions via DT, starting out from the RS so- 
lutions, we should use (|19|l combined with addition con- 
ditions (1201) and if^ . 

Now, let's take p9|) as a prop function. After the DT's 
(fTT)l execution we get 



where 



„(i) 



(1) 



v^'^'liy > 0) and v'l^ = v'^^^y < 0). This 



while V 

potential suffers a jump on the brane: 

.i^'(0)-.«(0) = -4M/3f# + °^+"^^ 



2n a,^ 



which is equally zero if the tension is 



(23) 



This case is really special for it allows to get the potential 



(1) 



from w^^^ by simple permutation of a^^ and 6^^ . Let's 
choose = l/-0i± [13 ■ It is clear, that for y — > 0: 

which, according to H15() is exactly the way it had to be. 
Using the tp^^ we receive a new metric 



It's asymptotical behavior is: A'^^\y) e ^^/^/y/o", 
when y +00, and A^^\y) — > e^'^^^/ ^b-^_ , when y 
—00. 123|1 is a particular case of a^^ = thus 



1 



±00. 



If we additionally insert condition = ^^i- = 0, we will 
get, after the substitution (3 —(3, the reduction of A'^^^ 
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to the solution of the RS set © , with the potential being 
equivalent to the initial V — — 3/i^/8. 
Returning to the general case, we receive 



2 ' 



and, after the simple calculation, we come to well-known 
sine-Gordon potential: 



3^2 



V3 



(0±(2/) - 0o±) 



where 



4'{y)± = 0o± + \/3arctan (^]J^ 



and 0o± are the arbitrary constants. 

Finally, to make the picture complete, we'll derive the 

quantity ("'o^^) ■ 

/z-\/3sin (a^ — «_) cos (a^ + q;_) . 



Note, that in the case H23fl we have 

(3 — 2/xcos2a+, a. 
where is a whole number. 



a+ + -{2N+l), 



IV. FURTHER DEVELOPMENT OF THE 
METHOD 



This paragraph is dedicated to three questions, each of 
them being connected with the further development and 
generalization of the method; namely: the n = 2 dressing 
of the RS solutions, the role of the shape-invariance and 
the two interesting ways of generalization of the model 
with the odd superpotential. The first way is realizable 
via the Adler theorem whereas for the second one the 
dressing chains of discrete symmetries should be used. 



A. Double dressing of the RS model 

Let's return to the RS model. The solution of it's 
Schrodinger equation with eigenvalue \2 — fi^ ~ can 
be chosen in form: 



V'.i = a.ie-y + K^e--y, ^Mv> 0), 

V'2- = V'i(y < 0). 



(24) 



Let's assume 

'^2- "^2+ + 2^ ' 

Using ini) for we get 



2v 



(25) 



^(1) ^ 

"± a^^e^y + h^^e-^^y 

r!± - (i. - m) (a,±a,^e(^+-)^ - h,J>,^e-^^^+-^y) + (^^ + m) {a,Ji,^e^-i^)y - a,^b,^e(^^-)y) 

I 



(26) 



Function H26|l will be used as a prop function for second 
DT. For the construction of a new ground state with 
energy A2 < we should choose > ^ > 0, a^^ > 
and < - this way the function -0^^' will not have 
the unwanted zeros. Consideration of H25I) gives us the 
conclusion, that value of positive coefficient a^^ should 
satisfy the inequalities: 



> 



2zy 



2iy + P 



/3>0 



02 

a. 



> 



2v 



2iy- P 



/3<0 



If these relations are correct, we are free to use (|16ll and 
()17|1 for the second DT, which results in: 



„(2) 



x2~2-f^logr!±. 




It is quite easy to make sure, that what we get is a brane 

(2) 

with correct jump condition and the tension ctq = +ao- 
Also note, that potential w^^^ has two levels: Ai = and 



B. Shape-invariant potentials 

It appears, that the proper work of DT method is only 
provided for exactly solvable potentials. One of the effi- 
cient ways of obtaining of such potentials lies in a usage 
of so-called shape invariants . If an initial potential is 
a function of y and some free parameters Ui : v = v{y; ai), 

and after the DT one gets u'-^-' = v'^^^y; a[^^) then v is 
called the shape-invariant (SP) potential. SP-potentials 
are quite common in quantum mechanics, e.g. the har- 
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monic oscillator The major point here is that ex- 
actly solvable potentials from P| (for the models with- 
out cosmological expansion) are also the SP-potentials. 
We'll confine ourselves with considering the three exam- 
ples from the cited article, retaining the terminology. 

A. Exponential potential. Superpotential is W{(j3) = 
2be~'^, where = log(a — by). Therefore 

[by - ay 

which is a well known SP potential: w^^^ = const x v. 

B. Even superpotential. For the DFGK model we get 
(with A{0) = 1): 

logA(y) = -^+5(l-e-2^^). 

In this case 

v{y) = lgb{2a-3b)e-''y + 16g'b'e-^'y + ^. 
After the DT 

V — > v^^^ = V — 4:-—r log A, 
dy^ 

we get 

y^^Hy) - lgb{2a + 3b)e-^''y + Wg^^e'^'y + 

o 9 

Thus f can be obtained from v by substitution a —a 
and g — > ~g. It means that v{y) is an SP-potential. 

C. Odd superpotential. In this case we have 

log^(y) = -ay-by'^. 

A calculation yields 

v{y) = 4{2by + af - 46. 

It is nothing else, but the harmonic oscillator and, hence, 
an SP-potential: w^^^ = v + const. Let's consider this 
example as the fuse in further constructing of the new 
exact solutions via DT. 



C. The generalization of the odd superpotential: 
DT on the harmonic oscillator background 

As we have seen, the odd superpotential W{(j){y)) — 
2K&^0(y), with K = -1-1 [y > 0) or K—-1 (y < 0) results 
in 

(l^iy) = 00 + Kb'^y, 

(27) 

logA(y)=logAo+ ^°'f^"\ 

where Aq = A{0), (/jq = (/'(O). Using (HTjl we get the 
Schrodinger equation of the harmonic oscillator 

Its ground state has the wave function ipi (y; k) = 
^-2{i^b\j+4>of/3^ and energy Ai = 4b* /3. Since 
V;i(0;+1) = Vi(0;-1) and 

V'i(0;+l)-7^i(0;-l) = -/3e-2*o/3, 

the tension will be /3 = 86^0o/3. All other cigcnfunctions 
and eigenvalues can be obtained via the formulas 

V'„(y;K) = (Sy -462(fe2y + «0j,)/3)"-i^^(y.^)^ 

(28) 

A„ = 462(2n- l)/3. 

In the article Adler has suggested the more general 
way of deletion of a groups of an excited states. It en- 
sures the general state A^,^^ of a transformed Hamiltonian 
to stay the same, though being described by a totally 
new wave function. Using this theorem and l|lt)|) one can 
delete the even number of adjacent levels (if the num- 
ber will be odd, for example N=l like before, then the 
new potential appears to be singular). Upon deletion of 
levels A3 and A4 one gets Higgs-like potential v^'^\y;K), 
with the ground state ^^\y; k) and the tension /J^^^: 



(2) _ 16b-^ 4096x^° + 12288a;^ + 21888a;'^ + 10368a;^ - 22599a;^ + 2187 (2) 8x^ + 3 _^^2^^ 

~ ~9 (64x4 ^ 27)2 ' ^1 ~ 64x4 -1-27*^ 

(29) 

(2) _ 8b^ 00 (51206 960^4 ^ 79202 „ 243) 
' ~ 3(6404 + 27)(802 + 3) 

I 



where x = nb^y + 0o. The potential v'^^^ is represented After the deletion of levels A3, A4, Ag and Ay we receive 
on the Figure 1. one very interesting potential which is shown on Figure 
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FIG. 1: The harmonic osciUator without levels A2 and A3. 
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FIG. 2: The harmonic oscillator without levels A3, A4, Ab, A7. 



2. Finally, the Figure 3 represents the potential which is 
obtained through the elimination of the levels A4 — A7. 
The expressions for the potentials eigenfunc- 
tions (including the ground state) and tensions can 
be easily obtained via the and H28|l. We omit them 
only for their extreme bulkiness. 

The reason why these potentials have such Higgs-like 
form is evident. In fact, the spectrum of potential with 
two symmetric minimums necessarily contains the lacuna 
between two down levels and all other spectrum. Delet- 
ing the levels A3 and A4 we construct the spectrum with 
such lacuna so it is only natural that this procedure re- 
sults in nothing else but Higgs-like potential. In a similar 
manner, elimination of the levels A4— A7 leads to the spec- 
trum with lacuna between the three bottom levels and 
the other spectrum. As a result one get the potential with 
three minimums (Fig. 3). Note that all these potentials 
have the following asymptotic: v^^^^y; n) 165^y^/9, 
when y — > ±00. 

Now lets return to potential l|29|l . Using the exact 

forms of the ti'^' and ■01^'' one can obtain the potential 
V{(j)) in parametric form (i.e. V — V{y) and (j) = (j){y))- 
The plot of the function d<j){x) /dx for y > is represented 
on the Fig. 4. This expression is real for | x |> ~ 
0.25735 so if we identify the </)o = then the field (piy) 
will be well determined all around the bulk. 
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FIG. 3: The harmonic oscillator without levels A4 — A7. 
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FIG. 4: The function 







D. Dressing chains and the fourth Painleve 
equation as the generalization of the odd 
superpotential 

Another interesting generalization of the odd superpo- 
tential (i.e. the harmonic oscillator) can be obtained with 
the aid of the dressing chains of the DT, which chains 
were previously introduced in l^^. Let's suppose, that 
we have defined the set of functions fk = dylogij/f!^ ^\ 

where ip^^ are n-times dressed solutions of the ^ . Then 
these functions will also be solutions of the dressing 
chain: 

ik + fk+i = fl - fi+i + au. (30) 

where are constants, and they can be expressed via 
spectral parameters Afc. Following we'll consider the 
periodic version fk = fk+N, otk = ctk+N with positive 
integer N which will be further referred as period. The 
whole theories of dressing chains are totally different for 
the odd {N = 2n + 1) and even {N = 2n) periods. The 
same is valid for the cases with a — and a 7^ 0, where 

N 

fc=i 

If a ^ and iV = 1 we get the harmonic oscillator. In the 
case A'' = 3 (a ^ 0) we get the fourth Painleve equation 
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(Pj-^^) for the function P{y) = — /i — y: 



P2 o 
P = — + -F3 

2P 2 



AyP^ + 2{y^-a)P+-, (31) 

where a = 2, a = — ai — a2/2 + 1, b = —0^/2. Thus 
the potential v{y) = /i + /i can be expressed in term 
of the Painleve-/y transcendents and has the foUowing 
behavior 



viy-yo) 



— -I- A/f7/cns I — ^ 

36 



+ My cos 1^-^ +(5o 



(32) 



Using H32|) one can conclude that this potential is indeed 
the generalization of the harmonic oscillator. In particu- 
lar, this potential is growing quadratically. 

We should also mention, that the spectrum of the v{y) 
can be found explicitly [231 . Moreover, it can be shown 
that the ground state is at Ai = 0. 



V. THE MODELS WITH COSMOLOGICAL 
EXPANSION 

If we are going to study the models with cosmolog- 
ical expansion, we should expose the DT method to 
some kind of generalization. This requirement inevitably 
follows from the fact, that the substitution = V 
now results in nonlinear equation rather then in linear 
Schrodinger equation. However, seemingly being as dif- 
ficult as it seems, this equation still can be reformulated 
in the way of linear spectral problem which does admit 
the Darboux transformation. In the next subsection we'll 
consider the spacetimes which are compact in y. We'll 
then introduce the simple orbifold 5*1 /Z2 and will show 
that, upon usage of DT, one can receive the solutions 
which result in exponentially small effective 4-D cosmo- 
logical constant. In the last subsection we'll briefly con- 
sider the new method which is valid for models both with 
and without cosmological expansion. 



A. Noncompact spacetimes 

If in 7^ then instead of ((HJ we have a much 

more complex system: 



2,^2 2ij2 2 



- i;^(^bS{y ~ yb), 



3 A 3 

6=1 

-^2 2.,_. 4P-2 



(33) 



3 3' ' A ' 
written for the set of N branes which are located dX y — 
yh- 

The function A" solves the following nonlinear equa- 
tion 



12H2 4 



A 3 



(34) 



6=1 



Therefore, we can formally introduce the spectral prob- 
lem (the Schrodinger equation): 



N 



Hy, A) = uiy) - f^bKy -yb)->^] ^{y; A), (35) 



fc=i 



where 



Uiy) 



;V{<i>{y);\) 



12i?2 



(36) 



We stress here, that U{y) is not a function of the spectral 
parameter A whereas V and -0 are. With this induced, 
we are free to use the spectral theory for the equation 
(p?5|l . For any given U{y) one can solve (p?5| to get ip 
and to find V from (|36|l . Of course, at the end of our 
solution we should, in order to find A = y^, receive i^j 
of positive value, but this problem is also avoidable, as 
has been demonstrated in sections 3 and 4. Indeed, one 
can construct positive solution via the DT starting out 
from any simple solvable model - taking the RS model for 
example. There is, however, one interesting difference 
with the case of stationary brane: the "kinetic term" 
appears to have the form 




-U + \\ - 



(37) 



therefore we can't just use the ground state solution like 
we did in the case of stationary brane, or at y — > oo 
wave function -0^0 and we will get (j)^ < 0. Obviously, 
this problem requires a more through examination. So, 
let's again start out from an "RS potential" U = /i^, 
and consider the case A = 0. The solution will have the 
form ifT^I) . Using it as the prop function in lfTT|l . one gets 

the same metric A^^^ from section 3, only with another 
potential. The main problem is now focused in the kinetic 
energy of the scalar field, having the form 



3H 



'^a^et^y + b^e^f^y. 



It is clear that for large enough y, one gets <; g. On 
the other hand, if (see l|2U) 



< 



sin 2a± 
~r ' 



then on the brane and in its vicinity, there are no such 
troubles. Therefore the simplest way of avoiding the 
problem is just dealing with the finite volume, which is 
less then the minimum value of y, resulting in the nega- 
tivity of the kinetic term. 

If we nevertheless wish to deal with infinite volume, we 
had to be sure to use only those positive solutions t^iy) = 
{ip+(y), at y > 0;ip-(y), at y < 0} which are growing 
as I y |— !■ cx). The construction of such solutions can be 
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done in different manners. For example, one can take 
the solution of the Schrodinger equation with A = A_i < 
A^^^ . The more interesting opportunity lies in applying of 
so called B-potential (Bargmann potentials) U{y) whose 
solutions of the Schrodinger equation all have the form: 
-0 — F{y) exp{ay) where F{y) is some polynomial |2.l| |. 
To show them in work lets rewrite the equation (|37|l in 
the form 



1 



1 



(38) 



where 



ip^/{2il}j_). The necessary criterion here is 



that ii^ < for any given y. Let's choose the in 
typical " Bargmann form" : 



= {Ay'^±By + C) e^^y , 



(39) 



where a, A, B and C are positive constants, associated 
by the condition 



2AC < < 4:AC. 



(40) 



The jump condition (|13|) . imposed to single brane at y ~ 
results in /3 = ~2{B + aC)/C . The solution of the 
equation H33(l appears to be 



A^ = e^''y/^^Ay^±By + C. 



(41) 



where A{y) = at ?; > and A{y) — A_ at y < 0. 
The condition H40|l guarantees that this solution will be 
nonsingular on brane as well as in bulk (since A(jj) ^ 
for any values of y) and also that ii^ < for y > and 
u_ < for ?/ < 0. Finally, choosing 



B^ - 2AC 



ACVC ' 
and using 1411) and 1)38(1 one get 

</>i _ 2A2y2 ± 2ABy + B^ ~ 2AC 
T ~ 4{Ay^ ±By + Cy 



^Ay^±By + & 



where 0(y) = (/)^ at ?/ > and (^(y) = at ?; < 0. 
Now we can see, that it is always possible to find large 
enough values of a that will greatly decrease the second 
member of equation with regard to the first, thus leaving 
the whole equation positive for both y > and y < 0. 



B. Orbifold model 

In previous sections, our efforts were wholly concen- 
trated on spacetimes that are noncompact in y. Now 
we'll consider the compactified case. As we shall see, the 
DT in this case allows to obtain the solutions, resulting in 
exponentially small 4-D effective cosmological constant. 

Let's choose the potential in l|35|l as V 
We'll start with the two solution of this equation: 



p? sin^ I 



■01 = cosh(/isin0j/), = sinh(/iy ± 6), 



where h =const. The function tl)^ is the solution of the 
(|35|l with Ai = whereas ?/'2± with A2 = cos^ d. 
Using as the prop function we can dress ■02± with the 
help of the DT The resuh wiU be 

'^tl ^ cosh(/iy ± 6) — sin Q tanh(/i sin By) sinh(/ij/ ± h). 

(42) 

where we should choose the sign " +" for y > and sign 
" for y < 0. One can show that -0^^^ > for all y. In 

this case the dressed potential U'^^^ will have the single 
bound state with A = Ai and eigenfunction 1 /'0i . 

Now we can compactify the y direction as an Z2 orb- 
ifold, with two branes sitting at each of the two fixed 
points (y = and y = L). The new aspect here is that 
S1/Z2 orbifold model permits only one bulk space. Due 
to the symmetry of the model, we only have to consider 
the jump conditions at y = and y = L. The direct cal- 
culations result in tension ctq = — 3/itanhfecos2 0/2 for 
the brane which is located at y = (we'll choose it as 
the visible brane) and 



3/1 sinh(/iL + b) (cosh^ (p sin 9L) — sin^ 9) — sin 9 sinh(/i sin 9L) cosh(/t sin 9L) cosh(/LtL + b) 
2 cosh(/isin0i) (sin^?sinh(/isin6'L) sinh(/LtL -I- 6) — cosh(/isin0i) cosh(/xL -I- 6)) 

I 



as a tensions on y ~ L branes. These expressions can be ^ / ^J~Mfi) ^ H/ ( ip^^'^ (0)) then 
ereatlv simolificd bv the choice b = —uL: V ^± / 



greatly simplified by the choice b = —fj,L 
(To = tanh(/iL) cos'^ 9, 



3/i 



tanh(/i sin 9L) sin 9. 



H{0) 



H 



^y/cosh{fiL) 



Since the expansion rate on the visible brane is H{Q) — Therefore, for large values of pL the effective 4-D cosmo- 
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logical constant A^f f on the visible brane is 



we get (for the case of single brane located at i = 0) 



A, 



ij2(0) 



-y/cosh(/iL) 



(43) 



This means, that 4-dimensional cosmological constant, 
as seen by observers from the visible brane, becomes 
exponentially small while the L grows large, therefore, 
automatically solving the cosmological constant problem 
right in the framework of the the brane worlds models 
(this idea was first suggested in article of S.-H.H. Tye 
and I. Wasserman Q). 

All we need now is to make sure, that (p^ > at 
< y < L (since we have the single bulk between two 
branes then it is enough to consider the case < y < L 
solely). Let's choose 6 ~ tt/2. In this case A2 ^ Ai = 0. 
We stress, however, that 6* 7^ 7r/2 exactly so A2 = —ii^a^ 
where a = 7r/2 - 6* < 1. To obtain 0^ (q) > it IS nec- 
essary that cosh(/iL) > IQH'^/ whereas (/)^ [L) > if 
cosh(/xL) < fi/{2H). The combination of these condi- 
tions is: 

— j <cosh(^L)<— , 

which is possible if /i > 2H . In fact, it is better to imply 
an enhanced condition /i ^ 2iJ, because large values of 
/i allows to avoid the conclusion that do = ifio? /2 ^ 0. 
Since /i is not bounded from above, we draw the wanted 
conclusion, that our model indeed permits positive (jP^ 
a,t < y < L. 

Nota Bene 3. Of course, the described theory 
could not generally protect some finite number of 
regions (at < y < L) from special situations when 
4'\{y) < (if ttiG value of ji is not sufficiently large, for 
example). In this case, however, one can use the method 
which was supposed in j24j . viz. these regions can be cut 
out, with their boundaries being sewn together in such a 
way that neither the warp factor (along with its first two 
derivatives) nor Ricci scalar will experience a jump. This 
unexpected fortune raises from the fact that matching 
is done at an inflection point of logyl (more correctly, 
this is approximately true if ij) is sufficiently large to 
neglect the term H/^/ijj in cither one of equations 
and ijSHIl)- 



Green function method 



-A + ( U - -ao(5(y) ) A = ~2H' 



(44) 



Solving linear nonhomogeneous equation 144|) for given 
U we can find 4>{y) and V{(j}) via the system 




Now, instead of the (|44|l let's consider the nonhomoge- 
neous spectral problem 



-7^(y; A) + ( U(y) - A - |ao<5(y) ) V(y; A) = -2H^ 



which can be written as 

L(y; A)V'(y; A) = -2i/2, 



(45) 



with L(y; A) = dP/dy-" + U{y) - A - 2n^6{y)/i. 

To solve H45|l we introduce the homogeneous linear 
equation 



L(2;;A)V'o(z/;A) = 0, 
and the Green function G{y, z;X): 

L(y;A)G(y,z;A)=%-z). 



(46) 



By solving these equation for given U{y) and A we can 
construct the solution of l|45|l as 

V'(y; A) = My, A) - J dz G{y, z- A). 

If i!{y\ A) will be positive then we will be able to define 
A — \pt\) just like before. But first we have to determine 
those conditions which make such positive solutions pos- 
sible. This is still the open question which will be con- 
sidered in the next publications. We should only note, 
that Green function can be constructed via solutions of 
the (|46|l . If initial potential U(j/) is exactly solvable one, 
then the same is true for N-times Darboux dressed poten- 
tial U^^-*. Thus, if we know the Green function G(2/, z; A) 
of the initial equation then we know the Green function 
G^^^{y, z; A) for equation with the dressed potential. In 
other words, the DT is appliable here, at least formally. 



This is another way to construct solution for the mod- 
els with cosmological expansion. First of all, we can di- 
rectly transform the system (jHHIl into the linear equation 
rather than use substitution xjj = A^. Indeed, introducing 
the potential 



VI. CONCLUSION 

This approach can be easily generalized to the case of 
if-branes. For this we should replace the term (35{y) in 
^ by 

^l3bS{y~yb)-, 

b 
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and also replace the two jump conditions H13() by 2K 
jump conditions at y = yi,. This is nothing more than 
technical matter, that's why we omit the details. 

Another one kettle of fish is that even in case of gen- 
eral position, where we cannot reconstruct the form of 
the potential V{(j)) as a function of 0, all these solutions 
will have the RS asymptotics (of course, if constructed 
precisely according to section 3). The DT gives a link be- 
tween the solvable problems, and probably most (if not 
all) of the exactly solvable potentials, from the poten- 
tial of the harmonic oscillator to the finite-gap potentials 
[25| . can all be obtained via these transformations. The 
physical sense of these potentials in the brane world is 
not clear. But our main purpose in this work was to 
merely demonstrate and advertise the DT as being the 
very powerful tool for manufacturing the exactly solvable 
potentials in the 5-D gravity with the bulk scalar field. 
It has been demonstrated that this method can be useful 
both in models with and without the cosmological expan- 
sion, although in the first case the situation is more com- 
plex and deserves the further studying. Notably, there 
is one imperfection of this method: the DT would work 
for the 5-D gravity only. If Z3 > 5 then the warp factor 
is a multi-variable function. Unfortunately, there is no 
cogent and effective general DT theory in several dimen- 
sions But, nevertheless, for Z? = 5 the DT is by far 
seems to be the best way to construct the nonsingular 
exact solutions. 
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